Functionally independent conservations laws in a quantum integrable model 
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We study a recently proposed quantum integrable model defined on a lattice with N sites, with 
Fermions or Bosons populating each site, as a close relative of the well known spin i Gaudin model. 
This model has 2N arbitrary parameters, a linear dependence on an interaction type parameter x, 
and can be solved exactly. It has iV known constants of motion that are linear in x. We display 
further constants of motion with higher Fermion content that are are linearly independent of the 
known conservation laws. Our main result is that despite the existence of the higher conservation 
laws, the model has only N functionally independent conservation laws. Therefore we propose that 
N can be viewed as the number of degrees of freedom, in parallel to the classical definition of 
integrability. 
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Introduction Quantum integrable systems have 
emerged from an esoteric beginning in the pioneering 
works of Bethe, Onsager, Yang, Baxter Ref. |l|) and 
others, and are of interest to a much wider commu- 
nity in the recent years. Optical lattices where quantum 
quenches can be realized to study dynamics away from 
equilibrium Ref. (0) Ref. |3), Ref. Q), transport theory 
Ref. (Q) Ref. (@) and considerations of entanglement en- 
tropy Ref. (0) are some areas where the the standard 
models of quantum integrability, such as the Ising model 
in a transverse field, the anisotropic Hciscnbcrg or XXZ 
model, and the 1-d Hubbard model find wide applica- 
tions, unanticipated in the pioneering studies. 

In this renewed era of interest, some basic ques- 
tions about quantum integrability remained unanswered. 
While the Yang-Baxter equations provide a deep mathe- 
matical underpinning to this field and indeed a source of 
most of the quantum integrable models, there are good 
reasons to look more broadly at the field. For instance 
there are other families of quantum integrable models 
that do not not fit naturally into this scheme, such as 
the Calogero Sutherland systems. Also, there is a well 
evolved language of classical integrable systems Ref. (@), 
with the dictionary and numerology of "degrees of free- 
dom" and of matching "functionally independent conser- 
vation laws" , that is not easily translated to the quantum 
arena. The meaning of the phrase degrees of freedom, 
is not always clear in many models of current interest. 
The number of conservations laws of even well estab- 
lished models, such as the XXZ model, is not quite a set- 
tled story, since new conservation laws have been found 
Ref. (0) in very recent studies, with important implica- 
tions for the transport behavior of this model Ref. (§)• 
Discussions of the meaning of quantum integrability is 
therefore of interest and has been addressed in recent 
literature Ref. Q), Ref. Ref. (dH) and Ref. (Q3). 

In this work, we study a recently proposed quantum 
lattice model Ref. (@) that is defined on a lattice of N 
sites, with either Fermions or Bosons, at each site, that 



is similar to the well known Gaudin spin | modelfil [l3|. 
For the Fermion model, we show that N, the number of 
lattice sites, can be viewed as the number of degrees of 
freedom in close analogy to the classical definition. To do 
this, we are able to enumerate the basic N conservation 
laws from the 1 particle sector. These basic conservation 
laws are linear in a parameter x, this linearity (or simple 
polynomial dep endence) has been highlighted in recent 
works Ref. (|12T ). Our most important result is as follows: 
we find that while other non trivial (i.e. four Fermi and 
six Fermi) conservation laws do exist in higher particle 
sector, even linearly independent from the basic ones, 
these are functionally dependent on the basic ones. We 
display the generating functions for the higher conserva- 
tion laws, and the non trivial functional relation between 
these. In the case of Bosons and the Gaudin model, using 
an interesting algorithm, we have checked that up to at 
least N = 20, there are no other non trivial conservation 
laws that are polynomial in the parameter x, and hence 
the N conservation laws available from the single particle 
sector are again the entire set. 

To summarize the models studies here let us note that 
a prototypical quantum integrable model can be built 
up in two steps: (Step-I) We identify a family of real 
symmetric commuting matrices H{x) and H(x) in N di- 
mension, that depend linearly on a parameter x: 



\H(x), H(x)) = 0, for all values of x. 



(1) 



This problem leads to several classes of solutions _that 
differ in the number of independent Hamiltonians H(x) 
that can be found for a given H(x). These are referred to 
as the Type-M with various values for M as in Ref. (fl3 ) 
Ref. (|15[). Type-1 represents a maximal set of such ma- 
trices with at most N commuting partners (including 
the identity matrix), where all such matrices can be ex- 
pressed in terms of the projection operators nij = 



2 



and 3-/V arbitrary real constants {ji,£i,di} as: 
H(x) = d i Z i . 



(2) 



2j = 7r« + a; > .(3) 
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The prime indicates the exclusion of the summed index 
with the fixed index. The commutation of the N basis 
operators [Zi, Zj] = is at the heart of this construction, 
and provide the N "higher constants of motion" or dy- 
namical symmetries, which depend parametrically on x. 
This is the maximal set of commuting operators, whose 
number equals the dimension of the defining space- N. 

(Step-II). The next step is to embed the operators 
of Eq. ((3]) into Fock space, thus Zi — > Zi by writing 
these in terms of quantum field operators acting upon 
suitable spaces. This procedure thereby gives rise to a 
quantum integrable models with Hamiltonians H —> H = 
J2i^i Zi. We may use spins (i.e. hard core bosons), 
canonical Bosons or Fermions to get models that coincide 
in the single particle subspace, but in higher subspaces 
are essentially different models. Let us note the Fermi 
model: 



7r7*( r a « + a\a r ) - "ff.n s - ^ 2 s n r 



, (4) 



where {a,, at} 
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6ij, n r = a£a r , and [Z r , Z s ] = for all 
1,2 N and for all values of x. Different com- 



muting operators, linearly dependent on the parameter 
x, are expressible as in Eq. ([2]). With spins in place of 
Fermions, this model reduces to the well known Gaudin 
model, generalized to include the coefficients 7$, and with 
(canonical) Bosons we get yet another model. 

In every higher particle number sectors, these opera- 
tors for any choice of the embedding Fock space have ob- 
vious descendent representatives that commute with each 
other and with H . For instance in the 2 particle sector 
for Fermions or the spin-^ Gaudin model, the Hilbert 
space has dimension of (^), where the above operators 
provide N commuting operators. These arc of limited 
interest to us. Our aim is to seek out operators that 
are linearly independent of these descendent operators. 
In this larger space we could for instance, independently 
construct type-1 matrices that are („) in number, - os- 
tensibly there's room for additional, as yet undiscovered 
conservation quantities that act in the two-particle sector 
but are null in the one-. 

In fact our main goal in this work is to examine the 
proposal that N is the maximum number of functionally 
independent commuting operators for H . Our search for 
functional independence starts with the more restrictive 
but technically feasible search for linear independence, 
followed by the test of functional dependence. 



H{x) = T + xl[T, S] + W 
H(x) = T + x([T,S] + W 



Algorithm for Finding Additional Conservation 
Laws: To uncover any such additional conservation laws, 
one can use an algorithm that we next describe. This al- 
gorithm can be implemented numerically for a relatively 
small size of the Hilbert space Af, as an illustration for 
Fermions TV ~ ( N ) where r is the number of particles. 
We used M < (g°) and in our studies of the Fermi, Bose 
and spin-i Gaudin models for , and the results are de- 
scribed below. 

The algorithm uses the matrix decomposition out- 
lined in Ref. (@) and Ref. (fT3) and the current- 
constructing algorithm outlined in jl5f . In Af dimensions, 
if [H{x), H(x)} = 0, it can be shown that there exists an 
antisymmetric matrix S such that 



(5) 



where T, T, W, W all mutually commute and, without 
loss of generality, we take them to be diagonal matrices. 
Given this decomposition Eq. ([T]) is satisfied if 



[T,S\,[T,S\ + W,[T,S] - W,[T,S] =0 (6) 



For a given H(x), S, T and W can be determined. Find- 
ing a conserved current H(u), i.e. finding a solution for 
the non-zero elements of diagonal matrices T and W, is 
then a matter of solving ("^J simultaneous linear equa- 
tions in 2Af unknowns, where the number of linearly in- 
dependent solutions is the number of independent con- 
servation laws [l6| . 

Fermions: The Four Fermi Conserved Current 
We applied the above algorithm to the matrices arising 
from the action of the Type 1 Fermionic Hamiltonian 
H = J^j djZj Eq. ([4]), in the two particle sector and 
found that there is a 2N — 3 member family of mutually 
commuting currents. N of these are just the action of the 
original currents Z r (see Eq. (j4])) in that sector. The N — 
3 other conserved quantities correspond to the particle 
action of four Fermi operators with a generating function 
of the form: 



Q(a) 
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2*-f{a-Ei){a-£j) 



(7) 



2 («- £ *)(a-£j)(a-£fc)' 



where <f>ijk = uJi 3 n k +n l nj'jl andw l:) = 7.477 (ajay+atdj)— 
n.ijj — njjf. Operators analogous to Eq. ([J} that are 
indexed by "r" can be found by taking residues 



= lim (a — e r ) Q(a), 



(8) 
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these are seen to be (i) linear in x (ii) are TV — 3 in number 
and (iii) linearly independent of Z r . The result (ii) re- 
quires an elaborate proof of three relationships between 
the apparently independent Q r . In particular one can 
show that Qi = 0, and we skip the remaining two for 
brevity. 

Analytical Proof of commutation: In what fol- 
lows we will prove that the generating function Q{a) is 
conserved i.e. 



[Z r ,Q(a)]=0 



(9) 



for all r = 1, 2, . . . , N and all values of x and a. Towards 
this end it is useful to define the quantity fiij = a\a,j — 

The commutator has terms of 0(x n ) with n = 0, 1,2 
and each must vanish identically. The O(x ) terms van- 
ish as all number operators mutually commute. 

To satisfy Eq. © to 0(x) we must show [I?} that 
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(a - Ei)(a - e k ) (a - e 3 )(a ~ s k ) 
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(10) 



(a - - Ej)(a - e k ) '' 



where we will subsequently follow the convention in which 
the prime on summation implies that all indicated indices 
are distinct. The RHS evaluates to 

_ y> iiijPijiik 

2-?{a- Ei)(a - £j)(a - e k ) 

In the LHS, we note that [wij,nin k ] = —[wij,rijn k ] = 
Pij n k and using the partial fraction identity 



(- 



Ei — Ej a — Ei a 



) 
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we see that Eq. (|10p is satisfied. 
To 0(x 2 ) we need to show 



E 



K-,Q'(c*)] 



0. 



(12) 



where Q'(a) = -£Q(a). We see from Eq. (0 that Q' 
depends upon three indices and so for the fixed pair in 
Wij in the Eq. (fT2"j). we may organize Q' as follows. The 
terms not involving either i or j commute and can be 
neglected. We may have both these indices in a class of 
terms Q', or only one index in another class of terms. 
Let us look at them separately: 

Two indices common: The relevant part of Q' in 
Eq. (TT2"!) may then be written as 



kjk + 4>ki 3 + 4>jki) 



E 



(a - Ei){a - Sj){a - e k ) 



Now Wij commutes with the first of these terms and for 
the others we calculate 

[wij,(t>uj\ - -ill + 7j) Pij n k + lilkPjkni - ijlkfiiknj, 

(13) 

and the term [wij,4>j k i] is obtained by exchanging if) j 
in Eq. (fT3"|) . giving the negative of the first so that the 
sum vanishes identically. 

One index common: The relevant terms in Eq. (fT2"]) 
may be written as 



Ei 

3 ,k<i 



4>ikl + 4>ilk + 4>kU 

/ V(" - - e k )(a - ei) 



(14) 

where we have used the symmetry of <j>ij k in the first two 
indices to isolate all terms with a triple of indices ilk that 
have one index i. It is easy to see that 

[wij,<Pku] = PijWki (15) 

[Wij,4>iki + folk] = -7i7j'7fe Wikjj - Pjkli) n-i- (k o I) 

(16) 

[Wij , 4>kli4>ikl + 4>ilk] = JiJjMijkt. (17) 

Here jijjMij k i is defined as the sum of the two terms 
Eq. ([T5|) and Eq. (flB)) . It is clearly symmetric in k -o- I 
and also antisymmetric in exchanging i ■(->• j, so that 
using Eq. (fTTj) we can rewrite the LHS of Eq. (fT2"]) as 



E 

j,k<i 



ji'-fjMi jk i 



(a - Ei)(a - £j)(a - £k)(a - £{) 



(18) 



We note an identity satisfied by cyclically permuting the 
indices jkl: 



JjMi jk i + j k M ik ij + jiMu jk = 0, 



(19) 



whereby the commutator vanishes identically. 

Functional relationship between the Generating 
Functionals: Let us note the generating function for the 
two Fermi currents: 



2 (a — Si) {a — Ej) 



so that [Z(a), Z(/3)] = for any a, (3, and 
Zi = lim (a — Ei)Z(a). 



, (20) 



(21) 



This is similar to the generating function for the four 
Fermi currents in Eq. ([7]). For completeness we note that 
from Eq. (0 



3,k 



{Ei - £j){£i - Ek) ' 



(22) 



To summarize, we computed a bilinear generating func- 
tion in Eq. (|20l) and showed that it has a four Fcrmionic 
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set of constants of motion Q r 's in Eq. <JSj> , with a gen- 
crating function Eq. (J7J, that are linear in x and also 
linearly independent of the Z r in Eq. ([2~Tj) . Remarkably 
these four-Fermi currents, while linearly independent of 
the Z's are not functionally independent of them. Indeed 
a brief calculation shows that 



Q(a) 



Z{af 



x{a) 




where x{a) = (£ m 7m/ ( a ■ 

To acquire some feeling for these generating func- 
tions, it is useful to transform to the diagonal represen- 
tation for the quadratic Fermion model. In Q one finds 
modes that diagonalize all the Z± and, by extension Z(a) 
by constructing the x-dependent canonical fermion crc- 
ation/annhilation operator set: 



4 = E 



Xi(x) 



with !'/,.</,! 



(24) 



CW( a ) = (1 - I ^ J ) ^™ } (a), then these satisfy the 
differential equation 



Though the commutativity of the generating functions 
is a straightforward consequence of Eq. (|27p. deriving 
the linearity in x of & m \a) is a bit more involved. 
In the (unrotated) original representation the linear x- 
dependence can be made explicit, i.e. 

C {m \a) = 



E 
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a - s k . 



2 = Ej (7j/(M«) - £ j)) 2 and 
Em^m/^it 1 ) — £ to) = V x - It follows that the two- 
Fermi generating function can expressed as 



where we require that 

v2 
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Z{a) 



1 



5; (a) 



E 



a - X m (x) 



(25) 



where n m = d^dm, x-dependent fermionic number oper- 
ator. Eq. (|23[) can be rewritten in terms of these opera- 
tors as: 



Q(a) = 1 - - 



x(a) 



E 



(a - Aj(x)) (a - Aj(x)) 



(26) 



In this representation, since both generating functions 
are constituted of the number operators in the diagonal 
basis, their mutual commutation is evident. However, 
known the linear dependence of the two generating func- 
tions on the parameter x is now hidden. 

Higher order Fermionic conservation laws: We 
can continue the logic of the above section to higher par- 
ticle sectors. While a more complete discussion is possi- 
ble, we shall content ourselves with a short presentation 
of our results. The Type 1 Fermionic model appears to 
possess a hierarchy of conservation laws, e.g. four-, six- 
,cight-Fcrmi terms, etc which are linear in x. We find 
from the above algorithm that a sequence of generating 
functions exist with m = 1,2,3,... where 



CM( a ) = i-(l 
m! 



x{a) 



e n 



n kj 



a - A fe . (x) ' 



(27) 

is a 2m-Fermi conserved current also linear in x, for 
m < N. Thus Z and Q discussed above correspond to 
m = 1,2. If we denote reduced generating functions by 



(29) 

where k\ , ki , . . . , k^ is an arbitrary permutation of the 
list of integers from 1 to N. For example, the six-fermi 
conserved current is: 
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ki,k 2 ,k 3 



1 n kl n k2 n k3 

6 {a - e kl ){a - e k2 )(a - e k3 ) 



l-x 
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s>3 



it 



E 



n kl n k2 (al 3 a ks +a^a fc3 
4 ^ (a- e kl )(a - s k . 2 )(a - e H )(a - e ks ) 



(30) 



The proof proceeds by induction. Assume Eq. (|29j) 
true for C^ m ~ l \a). Note that the product (1 - 
x/5(a))- 1 C( 1 )(o!) • C'( m - 1 )(a) contains (m - 1)-Fermi 
terms, as well as m-Fermi - we are concerned with the 
latter, as the former are exactly cancelled in Eq. (|28|l by 
the ■j^C ijn ~ 1 \a) term. The m-Fermi terms are given by 
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x(a) 
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x(a) 



mi a - e fc . 



(31) 
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(1) (m-l) (1) (m-l) 



m 1 



i(o;) 



2*£ 
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2(m- 1) 



(1) (m-l) 
c 2 L 2 



h± ,fc2 j - - - 



t I t m-l 

(" - £ fc, ) (" - £ fcm) a - £ fci 

(32) 



model considered. It appears the be the first case where 
one can establish the functional dependence between a 
"higher set" of conservation laws on a more fundamen- 
tal set, surprisingly both being linear in the parameter 
x. With this, we have a close parallel to the classical 
idea of integrability, and a well defined enumeration of 
the degrees of freedom that have so far been missing. 
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t it m— 1 

a k a k m + a k a k s TT n k 
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(a-£fcj(a -£fc m ) AA 



(33) 



where d m = (2(m — 1)!) . Summing these terms yields 
a (l-x/5(a))- 1 C( 1 )(a)-C'( m - 1 )(a) linear in x and, thus, 
also linear in x. 

Bosonic and spin models. We have used the above 
algorithm to study the higher particle number sectors 
for Bosons on a lattice with N < 6 and the number of 
particles 2 < n < 6 and find that there are no linear 
conservation laws beyond the Zj that are linear in x. For 
the spin-i Gaudin model, we have studied higher particle 
number sectors with arbitrary 7.;, where the operators 
analogous to Zj do exist. In that case we find that again 
there are no other conservation laws that are linear in x. 

Discussion The Type 1 fermionic model appears to be 
unique in possessing these 'extra' conservation laws linear 
in the parameter x. Not only do the Type 1 bosonic and 
Gaudin models not have them, numerical experiments 
with randomly generated Type M — M > 1 - fermionic 
and bosonic models similarly exhibit no such additional 
conservation laws. As we have shown, however, these ex- 
tra Type 1 conservation laws are functionally dependent 
on the N 'original' ones - the Zj. 

In conclusion, we have provided a non trivial exam- 
ple of the enumeration of dynamical symmetries in the 
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nontrivial results for an integrable matrix H(x). 
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